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ABSTRACT 


Differential equations that arise in astrodynamics are examined 
from the standpoint of Lie group theory. A summary of the Lie 
method is given for first degree differential equations. The Kepler 
problem in Hamiltonian form is treated by this method. Extension of 
the Lie method to optimal trajectories is outlined. 
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Introduction 


The study of problems in astrodynamics , as well as applied 
mathematics in general, leads to differential equations. When the 
dynamic of a problem leads to non-linear differential equations 
there is usually no systematic method of attack for their solution. 
The use of continuous transformation groups, in particular Lie 
groups, provides a general method to approach to global 
integration. 

A study was made of some elementrary problems in astrodynamics : 
the one and three dimensional Kepler problems. These problems are 
important because they have a formulation by the variational 
methods of mechanics and the resulting invariants of the Lie groups 
are constants of the motion. It is thus tempting to use Lie theory 
to study optimal control theory in astrodynamics since the this 
variational theory may be cast into a form that is an analog of 
Hamiltonian theory in particle dynamics. 


TheoryrThe Lie method of extended groups 


The practical application of Lie groups to the solution of 
differential equations(Bluman and 0016,1974) stems from the fact 
* that if the differential equation 


dx 

remains invariant under certain transformations, given by 
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■ ( then the equations can be arranged in such a way that the 

independent variables do not occur explicitly. The re-arrangement 

\ 

: involves a change of variable which are well defined by the 

invarance property. That is the invariance implies a change of 
variable 

T-C = tux, tj } • V ~ V<X, u) 




such that the differential equation becomes 

dii = ; § Cv.) 

d V J 

The kind of transformation required defines a continuous 
one-parameter group. That is a mathematical entity with the 
propert ies : ( 1 ) the the product of two elements is uniquely defined 
and is an element of the group, (2) the associative law holds for 
elements of the group, (3) there is a unit element, (4) and there is 
an inverse element. 

The basic proceedure for solving ordinary differintial equations 
starts from the infinitesmal parametric transformation 

■*, * * f C M ) 

and 

9.-- V e 

where G is a parameter and £ and ^ are the generators of the 
transformation. Demanding invariance of the differential equation 
and the integral invariants (that is constants of motion) and V 
leads to the operational equations 
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where ^ is the Lie operator defined by 
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and is its extension 

i 17 I^X + ' 1 ?^ + 

The generators are determined from (1) and the invariants u and v 
(in dynamics called canonical coordinates are calculated from the 
first order partial differential equations given by (2) and (3). 

The calculation process can be summerized as, given the 
differential equation / 1 apply the Lie operators 7L 'and ]£ to get 

r’ji =° 

linear partial differential equations in the generators. Use 
J U - o ^ “ 1. 

and solve the resulting characterstic equations to get the 
invariants u and v. Form the equation 

- (£) i'u.') 

av — 

gaining a reduction in order and a final quadruture if possible. 
Kepler Motion 

In the dynamics problems that follow focus is on the Hamiltonian 
formulation of the euqations of motion (Leach, 1981 ) , that is 
consider the first order system 

d. . ■aij dje - - ?— 

3e " d* 

where p and q are the canonical momenta and coordinates and H is 

10-5 







' t 


the Hamiltonian. 

As an example of a Lie group solution to a differential equation 
consider the one dimensional Hamiltonian defined by 

H = 

where p is the momentum, x the position of a particle and/ 4 - the 
gravitational constant. The equations cf motion are given by 

JL, = x - -p =o 


and 


JL t = V 




- O 


Impose invariance under one-parameter group transformations, that 
is under the extended Lie operator 



iX 

= o =o 



only 

surviving 

derivitives are 
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These values lead to the equations for the generators 

1 - V + -o 


(*) 


+ y ^ " ' V ** If + £ fs 

Solving equation(4) for ^ and putting into(5)leads to solutions for 
the generators given by 

£ - c*>**l*,,-f — O $ ® 

The integral surfaces of the partial differential implied by 

are generated by the integral curves of the equations 

- 6x - dfy _ dju - dv - djy 

i ~ o ~ -rr ~ o ~ o i 
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there are the solution-: 

ZL - A V » -p 

use of these coordinates and the equations of motion gives a 
relationship ’ ' ween u and v 

2 i 


- <? t-w.y ) = ~ VU / > 

d‘> ~ 


which integrates 
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or that 


A 


- sy* =■ ^ 


which means that E the energy is an invariant under translations in 
time. 

From 

Jtc 

v 


there results 


A 


d to 


= ctf- 


which can be integrated to 


t .<„ . -m 




a known result for one-dimensional Kepler motion. 

It is also possible to outline the three dimensional Kepler 
problem. The Hamiltonian will have the form 


H A 


K?l 
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The extended Lie operator 



where P and Q are vectoi 

a ‘" « + 7-% + t'-i + 

acting on the equations of motion 

JL, s a ° A, ' £- 

gives the vector partial differential equations 

-J + (i* - £ (£‘^1+^ t)~° 


and 



. Put y from (6) into (7) and solve for £ and ^ . The process leads 

to 3 groups: translations in time, rotations in space, and affine 
r translations in space and time. In turn these supply the 3 

invariants of motion, energy, vector angular momentum and the 
Laplace vector. 

Optimal Trajectories 

Consider the problem of extremizing a scalar performance index 

X fWl, 

subject to the constraints 

A ~ 

where u is a vector of control variables. This state equation 
conbined with boundry conditions defines the pioblem. For an 
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extremal the Hamiltonian defined by 


h_- z‘\f. 

must satisfy the necessary conditions 


•y _ c>fi 
‘ " 2»; 


i -- - i - !L 

> A l dV; 


(?) 


The Lie approach to (8) for the Kepler coasting flight car be 
outlined as follows: from the Hamiltonian 

A-5i 


H = (?«A i- /t 


R 

Writing the vector R, A as 

4>=a , , £» = A 


consider point transfromations in R and space 


R . = R+ *7, A ,= A+ ff + 

where <y are vectors. Apply the Lie and extended Lie Operators 
to find the generators and the invariants. 


Conclusion 


It has been shown that the use of the method of the Lie theory of 
.tended groups in the context of Hamiltons equations of motion 
yield the constants of motion. A method had been presented for use 

\ of Lie groups in optima] control theory, where hopefully with 

* 

. further work the non-linear equations will lead to invariants of 
f 

* the problem. 
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